home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b131.asm < prev    next >
Assembly Source File  |  1992-04-28  |  3KB  |  91 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.31    Unsigned Integer Divide
  8. ;The unsigned integer divide operation divides two 32 bit unsigned  integers.  The following code di-
  9. ;vides d0/d2 with  the resulting quotient in d0 and the remainder in d1.  
  10. ; 
  11. ;              Unsigned 32 Bit Integer                       Program    ICycles
  12. ;                Division of d0 = d0/d2                      Words
  13.     eor    d1,d1              ; clear d1
  14.     do     #32,dloop          ;32 quotient bits            2        3 
  15.     rol    d0                 ;dividend bit out, q bit in  1        1 
  16.     rol    d1                 ;put in temp                 1        1 
  17.     cmp    d2,d1              ;check for q bit             1        1 
  18.     sub    d2,d1    ifcc      ;update if less              1        1 
  19. dloop 
  20.     rol    d0                 ;last q bit                  1        1 
  21.     not    d0                 ;complement q bits           1        1 
  22. ;                                                          ---      --- 
  23. ;                                                  Totals:  8       133 
  24.  
  25.  
  26. ;The final  remainder is not produced.  This program  may calculate only the number of quotient bits re-
  27. ;quired and has variable  execution time.  
  28.  
  29. ;              Unsigned 32 Bit Integer                       
  30. ;            Division of d0 = d0/d1, d0>=d1                  
  31.     cmp    d1,d0    d0.l,d2.m
  32.     eor    d0,d0    iflo
  33.     jlo    divdone        ; divisor > dividend
  34.     bfind    d0,d0    d3.l,d8.l
  35.     jmi    dive2big        ;dividend has 
  36.                 ;32 significant bits
  37.     bfind    d1,d2    d0.h,d0.l    ;find # of quotient bits
  38.     movei    #32,d3
  39.     move        d2.h,d2.l
  40.     sub    d0,d2    d2.m,d0.l
  41.     inc    d2    d2.l,d2.h    ;compute loop iteration count
  42.     sub    d2,d3
  43.     lsl    d2,d1    d3.l,d2.h    ;align divisor
  44.     do    d2.l,divloop_fast
  45.     cmp    d1,d0        ;perform test subtract
  46.     sub    d1,d0    ifhs    ;if no borrow, do subtract
  47.     rol    d0        ;mult remx2, save quo. bit (borrow)
  48. divloop_fast
  49.     not    d0    d8.l,d3.l    ;flip inverted quotient
  50.     lsl    d2,d0        ;clean off any remainder
  51.     lsr    d2,d0
  52.     jmp    divdone        ;done
  53. dive2big    eor    d2,d2
  54.     do    #32,divloop_slow        ;same algorithm as 1st routine
  55.     rol    d0
  56.     rol    d2
  57.     cmp    d1,d2
  58.     sub    d1,d2    ifhs
  59. divloop_slow     rol    d0
  60.     not    d0
  61.  
  62. divdone    end
  63.  
  64. ;The final quotient is not produced.  This program may calculate only the number of quotient bits re-
  65. ;quired and has variable  execution time.  
  66. ;              Unsigned 32 Bit Integer                       
  67. ;            Remainder of d0 = d0 rem d1, d0>=d1             
  68.     cmp    d1,d0    d0.l,d2.m
  69.     jlo    divdone        ;divisor > dividend
  70.     bfind    d0,d0    #0,d2.l
  71.     jmi    dive2big        ;dividend has 
  72.                 ;32 significant bits
  73.     bfind    d1,d2    d0.h,d0.l    ;find # of remainder bits
  74.     move        d2.h,d2.l
  75.     sub    d0,d2    d2.m,d0.l
  76.     inc    d2    d2.l,d2.h    ;compute loop count
  77.     lsl    d2,d1    d2.l,d2.h    ;align divisor
  78.     do    d2.l,remloop_fast
  79.     cmp    d1,d0        ;perform test subtract
  80.     sub    d1,d0    ifhs    ;if no borrow, perform subtract
  81.     rol    d0        ;adjust remainder
  82. remloop_fast     lsr    d2,d0        ;align remainder
  83.     jmp    remdone        ;done
  84. dive2big    do    #32,remloop_slow        ;same algorithm as 1st routine
  85.     rol    d0
  86.     rol    d2
  87.     cmp    d1,d2
  88.     sub    d1,d2    ifhs
  89. remloop_slow     tfr    d2,d0
  90. remdone    end
  91.